The local content of all pure two-qubit states 
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The (non-)local content in the sense of Elitzur, Popescu, and Rohrlich (EPR2) [Phys. Lett. A 
162, 25 (1992)] is a natural measure for the (non-) locality of quantum states. Its computation is in 
general difficult, even in low dimensions, and is one of the few open questions about pure two-qubit 
states. We present a complete solution to this long-lasting problem. 



John Bell showed in 1964 [l| that the "classical" meta- 
physical assunrptions used by Einstein, Podolsky and 
Rosen (EPR) 0, Q have empirical consequences in the 
form of a constraint on correlations of quantum mea- 
surement outcomes: the Bell inequalities. Bell later sub- 
sumed these assumptions into the notion of local causal- 
ity [1] , today commonly abbreviated to just locality. The 
experiments done so far give strong evidence that Bell in- 
equalities are violated as predicted by quantum mechan- 
ics and thus that there can not be any local model for 
the corresponding correlations: the world is in this sense 
nonlocal. 

Bell's theorem is a metatheoretical result of, at first 
sight, rather philosophical concern. That it can be of 
actual practical interest came as a surprise @. Starting 
with an intuition of Ekert in 1991 0, it was established 
later on that nonlocal correlations can be made the basis 
of so-called device independent quantum information pro- 
cessing There are many examples in the history of 
physics where a theoretical result is used for applications 
not foreseen by the pioneers, but Bell's theorem with its 
foundational character might be the most striking yet. 

Whereas the violation of a Bell inequality implies non- 
locality, it does not quantify it. To take the amount of 
the violation as a measure is problematic. E.g., the state 
which violates it maximally depends on the particular 
inequality fl3l |. A more natural measure of nonlocality, 
proposed three times independently jl4| - |T6i . quantifies it 
in terms of communication complexity [l7|: how many 
bits of communication are minimally needed for the sim- 
ulation of quantum states. For the simplest case of pure 
two-qubit states most answers are known. Surprisingly, 
as shown by Toner and Bacon [l^ , a single bit is already 
enough for the maximally entangled state and two bits 
for all partially entangled ones. The main open question 
is whether one bit is already enough even for partial en- 
tanglement or whether the two bits of [l^ are necessary. 
Whatever the answer to this question, communication 
complexity is clearly a very coarse measure for the non- 
locality of these states. In that respect a more promising 
approach starts with the so-called local content. 

Consider an experiment where one measures the proba- 
bility distribution Pq{A, B\a,b) of the measurement out- 
comes A and B for measurements a and b by prepar- 
ing repeatably particle pairs in the same quantum state. 
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FIG. 1. An EPR2 decomposition holds for all measurement 
directions a and b and measurement outcomes A and B. The 
red wavy lines indicate nonlocality. The nonlocal part is more 
nonlocal than the decomposed quantum probability distribu- 
tion. The outcomes of the local part are determined by A and 
the corresponding measurement direction. 



If Pq comes from a nonlocal state the whole ensem- 
ble of particle pairs cannot be described with a local 
model. However, as first observed by Elitzur, Popescu, 
and Rohrlich (EPR2) [ll, there are states for which it 
is possible to regard a fraction pL €]0, 1] of the parti- 
cle pairs as behaving locally, and the completing fraction 
1—pL as behaving nonlocally. We will then say that there 
exists an EPR2 decomposition into a local part Pl with 
weight PL and a nonlocal part Pnl with weight I — pl 
(see figure [1]): 



PQ=PLPL + il-pL)P. 



NL, 



(1) 



where the decomposition holds for all von Neumann mea- 
surements and all outcomes. EPR2 decompositions are 
in general not unique. What is unique is the local content 
PLC which is defined as the maximal value pl can take 
over all EPR2-decompositions. We will refer to a decom- 
position with PL = PLC a maximal EPR2 decomposi- 
tion and define the nonlocality of a state (or its nonlocal 
content) as 1 — plc- This nonlocality measure can be 
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defined for any probability distribution. For instance, 
the nonlocality of the nonlocal part of a maximal EPR2 
decomposition is 1 and such a nonlocal part is thus max- 
imally nonlocal. We also note that the nonlocal part of 
any EPR2 decomposition of a quantum state is always 
more nonlocal than or equally nonlocal as the quantum 
state it decomposes (sec figure [1} and is in general not 
quantum; but it is still non-signaling, since both Pq and 
Pl are. 

The local content of pure two-qubit states was already 
investigated by EPR2. By making appropriate choices of 
the basis, all these states can be written as 



(2) 



l*W> =cos( - 1 |00)+sin( - ) 111) 



with 6 £ [0, 7r/2]. All entangled states (6* > 0) are nonlo- 
cal [H [111. In spite of this, EPR2 found in 1992 an ex- 
plicit EPR2-decomposition with local weight — sin 6*) 
for all these states, thus giving a lower bound on their 
local content. Furthermore, they gave an argument that 
the maximally entangled state has a local content of zero 
and is thus fully nonlocal. Further progress was achieved 
only in 2006 by Barrett, Kent, and Pironio [l^. Since the 
local part necessarily obeys Bell inequalities, this places 
constraints on any EPR2 decomposition and can be used 
for deriving an upper bound on the local content. With 
this idea, and using a chained Bell inequality [1^, [l^] , 
they were able to rederive (without a technical assump- 
tion used by EPR2) the result that maximally entangled 
states are fully nonlocal [25|. Next, Scarani [la] narrowed 
the set of possible values oipLC further down: He derived 
both a new lower bound of 1 — sin0, and numerically, us- 
ing again the same chained Bell inequality as [24I, an 
upper bound of cosO; he conjectured in the preprint ver- 
sion [i^l of [2^ that this equals the local content, i.e. 
that PLC = cos 6*. This conjecture was proven to be cor- 
rect by Branciard, Gisin, and Scarani [28| for the states 
with cos 6* e [0,|]. 

The main result of this article is an explicit model for 
the local part which generalizes the model of [2^ and 
leads to a valid EPR2 decomposition with local weight 
Pl = cos 9 for all states ([2]). Since cos6' is an upper 
bound for the local content (which we prove analytically 
in Appendix |B|, this implies that plc = cos 9 and thus 
proves Scarani's conjecture for all pure two-qubit states. 

a. The problem We will parametrize the von Neu- 
mann measurements of Alice and Bob by two three- 
dimensional unit vectors a and b on the Bloch sphere, 
in the form [23 




^bj_ cos/3> 
b = I 6_L sin f3 



(3) 



With the definitions s := sm9, c := cos9, and x '■= ct — /?, 
the expectation values of Alice and Bob's measurement 



outcomes A,Be {+1, —1} for the states ([2]) are 

{AB)Qia_h = azbz + sa±b± cos X, (4) 

(^)Q|a = ca^, (5) 

(i?)Q|5 = c6,. (6) 



Without loss of generality |30|, |31| , we only consider de- 
terministic models for the local part. That is, we take the 
outcomes to be functions of some shared random variable 
A and of the corresponding measurement directions: 



A^A{a,X), B^B{b,X). 



(7) 



Because of locality, A is independent of 6, and B inde- 
pendent of a. Moreover, the probability density p{X) of 
the shared random variable is assumed to be independent 
of both measurement directions. The expectation values 
of the local part are then given as 

{AB)L\a.b = I dXp{X)A{a,X)B{b,X), 



(8) 



{A)L\a = J dXpiX)Aia,X), 
{B)L\b^ J dXp{X)B{b,X). 



Ultimately, our aim is to prove that p^c = c. Since 
PLC < c, as proven numerically in [2^ and analytically 
in Appendix iBl it is enough to find a model for the local 
part that admits a local weight pl = c. 

Barrett, Kent, and Pironio [24I showed that the max- 
imally entangled state, being fully nonlocal, necessarily 
has random marginals. Since the nonlocal part of a max- 
imal decomposition of the states ^ is also fully nonlocal, 
its marginals have to be random, too (we prove this in 
Appendix [C]) . As can be seen from ([5]) and ([6]), this im- 



plies the constraints (A) ^ 



and (B) 



6, for the 



local part. Note that this agrees with the ansatz used in 

Once one has found candidate functions ^ for the 
local part, it remains to be checked, whether the implied 
Vl = ^^^-r'^ (assuming = c [H]) is a 



nonlocal part Pj 



probability for all measurement directions and outcomes. 
Since the nonlocal part has random marginals, this is the 
case if the correlation term is between —1 and 1. The 
problem to be solved reduces to the following task [28| . 
for any given value of c: 

Problem. Find functions A{a, X) and B{b, A) taking 
values in {+1,-1}, and a probability distribution p(X) 
such that 

{A)L\a = a., (9) 

{B)L\b = (10) 

|(^S)Q|„,5-c(AB)i|„,b| < 1-c, (11) 

hold for all measurement directions a and b. 

Note that this formulation of the problem is redundant, 
since, as shown in Appendix [Cl constraint (|lip implies 
(HD and (UnD. 



3 



b. The model. Our star ting point is the model for 
the local part introduced in [2^, which we will call the 
BGS model. The complete model below will be a gener- 
alization of this model. 

BGS model for the local part. Alice and Bob 
share the random variable A which is a random three- 
dimensional unit vector uniformly distributed over the 
Bloch sphere. For measurement directions a and b, Alice 
and Bob output, respectively, 



A{a, A) = sgn[ar — a • A], 
B(6,A) =sgn[6, -b- A], 

where sgn is the sign function (Cz {+1, —1} ). 



(12) 
(13) 



This model gives the right marginals and ([T 
and the correlation becomes (see [23| or Eq. (jE15|) in 
Appendix [El for the explicit form of E^'^^ {az,bz, x))- 



(AB) 



BGS 
L\a,b 



( l-\az-bzl ifx = 0, 

<^ -l+|a,+6,|, ifx = ^, (14) 

I i?fG«(a„6„x), ifO<|x|<7r. 



This model solves the problem partially: as verified nu- 
merically in [2^, PT|) is fulfilled for all states such that 

ce[o,|]. 

To get a leading idea for our generalization of the BGS 
model, we will concentrate on a necessary condition for 
the local part [28|. If Alice measures in a direction a 
and gets the outcome +1, Bob's qubit collapses to the 
state |a(l)) (the eigenstate of a(l) • cr with eigenvalue 
-1-1), where \a(l)) is one step above \a) on the "Hardy 
ladder" [H, [SJ] (see Appendix The vector a(l) is 
given by: 



^a^{l) cos 
a(l) = I a±(l) sma | with az{l) 
az{l) 



1 



(15) 



and a±{l) = -^1 — az{l)'^. If Bob now measures in the 
direction a(l), he will always get the outcome -1-1 and 
thus Pq{+1, — l|a,a(l)) = 0. In a similar way, we also 
have Pq(— 1, +l|a, a(— 1)) = 0, where 



^a±{~l) cos Q;^ 
a(— 1) = I ax(— l)sinQ; | with az(— 1) 



caz 



(16) 



and a(— 1) is one step below a on the Hardy ladder. Since 
Pq = plPl + {^-Pl)Pnl, this implies (for pL ^ 0) that 
the corresponding probabilities for the local part must be 
zero, too: 



PL(+l,-l|a,a(l)) = 0, 
PL(-l,+l|a,a(-l)) = 0. 



(17) 
(18) 



The BGS model fulfills these conditions for all states, also 
for the ones which violate constraint ([TT|) and thus for 



which the model does not give a valid EPR2 decomposi- 



BGS 



tion. To see this, we note for instance that {A) 



i^)L\afi) — (^zi^), and that, since a and a(l) have the 
same azimuth angle (i.e., x = 0) and 0^(1) > a^. 



\az~az{l)\ = l + az-azil). (19) 



Hence, one indeed has 
Pf«^(+l,-l|a,a(l)) 



BGS 



(B) 



BGS 
L|a(l) 



(1) 



(20) 



Let us now see in what way we can alter Alice's output 
function in the BGS model, such that conditions 
(ITTI) and (fT8)) still hold (we ignore condition ([9|) for a 
moment). Assume that Bob measures in direction b = 
a(l): Condition ([T7| still holds if Alice chooses to output 
A{a{t),X) with 



^a^{t) cos 
a{t) = I ax(t) sina | , 

CLz{t) 



(21) 



and —1 < az{t) < az{l). Similarly, assume that Bob 
measures in direction b = a(— 1): Condition (|T8)) still 
holds if Alice chooses to output A{a{t),X) with az{—l) < 
Oz{t) < 1. Taken together, since Bob could be measuring 
in either direction a(l) or a(— 1) about which Alice has 
no information, Alice can choose to output A{a{t),X) 
with az{—l) < az{t) < az{l) and still be sure that con- 
ditions ([T7|) and (US]) hold. We will say that Alice is free 
to choose a{t) in the Hardy sector Ha of a. To use this 
wiggle room is precisely the idea for the construction of 
our 6'-dependent generalization of the BGS model. 

Bob's output will be identical to his output (fTS]) but 
Alice will now output A{a{t),X), where t e [—1,1] 
parametrizes the vectors in the Hardy sector Ha, such 
that a(0) = a, and that a(l) and a(— 1) are the two 
vectors introduced in ([TSl - ITB)) . Alice picks a particular 
a{t) depending on a local (nonshared) random variable 
t, which has a probability distribution Pa{t) that depends 
on Oz (besides depending on 6) [s^. One obvious con- 
straint on Pa (t) comes from the condition that ^ should 
still be fulfilled. Note that for x = 0, the model implies 
the same correlation term p4)) as the BGS model if b 
does not lie in the Hardy sector Ha of a. Incidentally, 
for such measurement directions, one can show (see Ap- 
pendix|D]) that the BGS model does lead to a vahd EPR2 
decomposition for all values of 9, as it should be for the 
following generalization having a chance to be valid for 
all states and all measurement directions. 

Complete model for the local part. Alice and Bob 
share the random three-dimensional vector \, uniformly 
distributed on the Bloch sphere. Alice can sample an ad- 
ditional local random variable t G [—1,1] according to 
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the following distribution, which depends on her setting 
a through its component : 



Pa{t) 



sln(7) 1 + G{t)a^ 



4c VI - G{t)^ ' 



(22) 



1 + c Y -I 

with 7:=T— , G(i):--T-T- (23) 
1 — c 7+1 

^/zce OTirf Bob 's outputs are then defined as 

A{a, A, <) = sgn[a,(t) - a(t) • A], (24) 
-B(6,A) =sgn[&,-b- A], (25) 



where 



a±{t) cos q} 



a, + G{t) 



a{t) := aj_{t)sina with az(t) := , ^, 

a,(t) / 1 + G(t)a, 



, (26) 



and a^(t) = ^Jl — az{tY 



The parametrization of the vectors in the Hardy sec- 
tors Ha as a{t) is motivated in Appendix [Xj The prob- 
abihty distribution Pa{t) is derived in Appendix lEl from 
the postulate that constraint ([TT|) is saturated for the case 
that b hes in the Hardy sector T-La of a. This, together 
with the previous remark on the case where b does not 
he in "Ha, aUows us to prove analyticahy that our local 
model above defines a valid EPR2 decomposition, with a 
maximal local weight = Plc = c, for all measurement 
directions a, b with x = 0. One can show in a similar 
way that this also works for measurements with x = tt. 

For measurement settings a, b with < |x| < tt, we 
checked numerically with standard numerical optimiza- 
tion tools of a common computational software program 
that constraint (jlip holds within the expected numerical 
precision (of absolute order 10~^°). Note that for any 9, 
(AB) j^^^ J, depends only on the three parameters a^, bz, 
and X- We also note that the implied nonlocal part is 
nonquantum for aU c € ]0, 1[ (see Appendix |F|. 

c. Discussion. Our new EPR2 decomposition for 
two-qubit states, inspired by that of [28| and based on 
properties of Hardy ladder vectors, thus allows us to show 
that the local content of all pure two-qubit states ([2]) 
equals c = cos 6. This proves a conjecture formulated 
in [i^l ■ Interestingly, the nonlocality 1 — c of these states 
is found to be monotonously related to its entanglement 
— against the trend of results which suggest that entan- 
glement and nonlocality are of a rather different nature 
(see e.g. (Mil). 



One obvious next step is to consider more general 
states. For the case of higher dimensional partially en- 
tangled two-partite states, we only know of the numerical 
lower bound on pLC found by Scarani [2^ for a class of 
entangled qutrits. We believe that the methods presented 
in appendices iBl and [Q should be generalizable to get an- 
alytical upper bounds which might then pave the way for 
models of the local part. Whether the monotonous con- 
nection between entanglement and nonlocality then still 
holds is a question of particular interest. 

We would like to end with one possible application 
of our decomposition and some ensuing open questions. 
Any EPR2 decomposition can be used as a starting point 
for a simulation protocol, where only the nonlocal part 
remains to be simulated [39|. Now, the arguably biggest 
difficulty in the simulation of partially entangled states 
is their nonrandom marginals. Using our maximal EPR2 
decomposition seems in that respect ideally suited, since 
its nonlocal part has random marginals. We also note 
that nonlocal correlations with random marginals that 
can be simulated with one bit of classical communication 
can be simulated with one nonlocal Popescu-Rohrlich 
(PR) box [iOl, too [U. Thus, if the nonlocal part of 
a maximal EPR2 decomposition for two-qubit correla- 
tions can be simulated with one bit of communication, 
then the corresponding quantum state can be simulated 
with one PR box. This is however known to be impossi- 
ble for weakly entangled 2-qubit states [s^; hence their 
nonlocal part in a maximal EPR2 decomposition can- 
not be simulated with one bit of communication. Does 
this answer the open question, whether partially entan- 
gled states can be simulate with a single bit? Not yet: 
Whereas it is always possible to simulate the decomposed 
state by simulating the nonlocal part, this might not be 
the most optimal way. 

On the one hand, maximal EPR2 decompositions lead 
to a natural alternative measure of nonlocality. On the 
other hand, they thus lead to new insights on the quan- 
tification of nonlocality in terms of communication com- 
plexity. 
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Appendix A: Hardy ladder vectors 

The vectors a(l) and a(— 1), one step above and 
below a ~ (a^ cos a, a_L sin a, a^) on the Hardy lad- 
der, are defined such that Pq(-|-1, — Ija, a(l)) = and 
Pq{—1, -t-l|a,a(— 1)) = respectively, and given by 

(aj^(l)cosa\ , 
a_L(l)sina with ajl) = -^^^^ -, (Al) 
a.(l) / ^"^ + 1 

and 

/a^(-l)cosa\ 
a( — 1) = aj^( — l)sinQ! with az{—l) ~ — . 

(A2) 

The successive steps of the Hardy ladder are iteratively 
defined as 

PQ{+l,-l\a{n),a{n+l)) ^0, and (A3) 
PQ(-l,+l|a(-n),a(-n-l)) =0 (A4) 



6 



for all n G N. To get the explicit form of a{n) for all 
n €E Z, we note that the maps 



a, ± c 



±ca. 



1 



(A5) 



are (real) Mobius transformations. The composition of 
Mobius transformations basically reduces to the product 
of their associated matrices, in our case to the product 
of 



1 ±c 
±c 1 



(A6) 



Defining (we assume c < 1) 



7 := 



1 + c 



>1, Gin):-- 



7" - 1 



-e[-l,l], (A7) 



1 - c ~ ' ' ' 7" + 
one gets the explicit form of a{n) as 

' a_L (71) cos a\ , N 

a[n) = a^(n)sma with az[n) := — — — , 

a,(n) / G(n)a, + 1 



(A8) 



and a±{n) ~ -^/l — az{n)^. This can independently be 
checked e.g. by induction. 

In order to get a parametrization of the vectors "be- 
tween" a(— 1) and a(l), wc let n take real values t in 
the interval [—1,1]. Note that a^(i) is a monotonically 
increasing function. We then define for every vector a its 
Hardy sector Ha to be the set of all vectors which can 
be written as a{t) for some t e [—1, 1]: 



Ha ■.= {a{t)\te [-1,1]}. 



(A9) 



Appendix B: Analytic proof that c is an upper 
bound for plc 

Let us label N measurement directions of Alice as 
a(-iV), a{-N + 2), . . . , a{N - 2) and N directions of 
Bob as 6(-7V + l), b(-iV + 3), b{N-l). For a given 
probability distribution P with binary outcomes A, B, 
we define as a correlation measure the expression 

In{P) =P{A-N B-N+l) + P{A-N+2 ^ B-N+l) + ... 

+ P{An-2 ^ Bn-i) + PiA.N = -Bjv-i), 

(Bl) 

where P{Ai 7^ Bj) is a shortcut notation for P{A 7^ 
B\a{i) , b{j)) . As shown in [2^, [^J], In is a chained 
Bell polynomial, defining the following chained Bell 
inequality — satisfied by all local correlations: 



iNiPL) > 1- 



(B2) 



Since by linearity. 



and since /7v(-P/vl) does not have any constraints except 
that the probabilities are between zero and one and thus 
In(Pml) > 0, In{Pq) is an upper bound for pl [22j : 



PL < In{Pq)- 



(B4) 



Wc show now that the measurement directions of Alice 
and Bob can be chosen such that 



lim In{Pq) = c. 

For that, we start with an arbitrary unit vector 

^v±{0) cos(a;) 
v{0) := ( wj^(0)sin(cj) 



(B5) 



(B6) 



for any ^^(O) e] - 1, 1[, v±{0) = yj\ - v^{^Y, and a; G 
[0,27r[. Let us then define, with the notation of (jA8[) . 



'v\_{n) cos(cj)^ 
v(ri) : = I v^(n)sin(cj) 



(B7) 



We then choose Alice's settings to be a(n) = v(ji) for 
n € {-iV, -N + 2, . . . , - 2}, and Bob's settings to be 
h[n) = v{n) for n e {-N + 1, -N + 3, . . . , N ~ 1}. An 
example of the z coordinates of such a set of vectors is 
drawn in figure [U for N = 3. Note that we have 



lim v(n) 



and lim v{n) 



(B8) 




b.i-2) 



FIG. 2. Hardy ladder for A'^ = 3. Every colored line represents 
one term in the chained Bell inequality, the dashed red one 
being the last one, P{A-m ~ Bn-i)- 

Since 



iNiPq) = PlIn{Pl) + {1~Pl)In{Pi 



NL 



(B3) 



PQ(l,-l|a(n),b(n+l)) = 



(B9) 
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for two neighboring Hardy steps a{n) = v{n) and b{n + 
1) = v{n + 1), one has 

{AnBn+l)Q (^-B)Q|a(n),b(n+l) 

= l + c[?;,(n)-w4n + l)], (BIO) 

and hence 



[v,in + I) - v,in)]. (Bll) 



Similarly, 



Pq(A„+i ^ S„) = + 1) - v,in)]. (B12) 

In{Pq) then becomes 

= %,{n+l)-v,{n)] + -[l + {A^NBN-i)Q] 



n=-N 
C 



^ MN - 1) - V,{-N)] + - [1 + 



And because of (IB8I). we have 



lim In{Pq) = c, 



which leads with (IB4I) to 



PL < c. 



(B13) 



(B14) 



(B15) 



Appendix C: Proof that (A)jy^ = (B)jy^ = for a 
maximal EPR2 decomposition of Pq 

As shown in Appendix [B] there exist measurement di- 
rections for which 



lim In{Pq) = c. 



(CI) 



Since we have, for a maximal EPR2 decomposition of Pq 
(i.e., for PL =Plc = c), 

iNiPq) = cIn{Pl) + (1 - c)lNiPNL), (C2) 

with In{Pl) > 1 and In{Pnl) > 0, then we must have 

lim In{Pl) = 1 and lim In{Pnl)^0 (C3) 

Af— >-oo N^oo 

for such a set of measurement directions. 

Note that apart from the last term P(A_jv = Pw-i), 
a term that is present in In{Pnl) is also present in all 
subsequent polynomials lN+2k{PNL), for /c S N. Since all 
terms appearing in In{Pnl) are probabilities, and since 
livaiq^oo In (Pnl) = 0, then all terms have to be zero, 



except possibly for the last one which must only tend to 
zero (as N ^ oo). 

Take one term: P{An ^ Bn±i) = implies that 

P(A„ = +l,S„±i = -l) 

= -(1 + (Ai)wL - {Bn±l)NL - {AnBn±l)NL) = 0, 
P(A„ = -l,P„±i =+1) 

= ^(1 — {An)NL + {Bn±l)NL — (A)iPri±l) JVi) = 0, 

and therefore 

P(A„ = +1, P„±i = -1) - P(A„ = -l,Bn±l = +1) 

1 



>NL 



{Bn±l) 



NL) 



0. 



(C4) 



Hence, {An)NL = {Bn±i)NL, and by induction {Ai)NL = 
{Bj ) for any odd value of i and even value of j . 
In a similar way, P(A_Ar = P^v-i) implies that 

P{A^N = Bn-i = +1) 

= ^(1 + {A^n)nl + {Bn-i)nl + {A_nBn_i)nl) ^ 0, 
P(A_jv = Bn-1 = -1) 

= ^(1 — {A-n)nl — {Bn-i)nl + {A-nBn-i)nl) — > 0, 
and therefore 

P{A^N = Bn-1 = +1) - P{A^N = Sat^i = -1) 

^li{A-N)NL + {BN-l)NL)^0. (C5) 

Now, from the previous argument, (A_jv)jvl = 
{Bn-i)nl = {A.,)nl = (Pj)iVL for any i,j. Eq. ([05]) 
then implies that (Ai) ml ~ {Bj) nl ~ for all the mea- 
surement directions appearing in In{NL). 

Since we started the construction in Appendix [B] with 
an arbitrary measurement vector (sec equation (|B6|) ). we 
must have, for the maximal EPR2 decomposition consid- 
ered here, 



NL\b 



= 



(C6) 



for all a and b. This in turn implies that the marginals 
of the local part are, as given in Eq. (PI ITU)) . 



{A)L\a = az, {B)L\b = b^. 



(C7) 



Appendix D: Proof that the BGS model gives a 
valid EPR2 decomposition for measurement settings 
such that X = 0, and b ^ Ha 

In this Appendix we prove that when Alice and Bob's 
settings are such that x = 0, but such that b does not 
lie in the Hardy sector Ha of a (i.e., > Ozfl) or < 
Q.z(— 1)), then the BGS EPR2 decomposition is valid 



for all states. For that, we will show that pT|) holds for 
all values of c. 



Since b is assumed to lie in Ha , bz can then be written 



Suppose indeed that x = £ind > 0^(1) 



> 



Qz- Then, omitting the conditional settings, we have 
{AB)q - c(AB)l = cizbz + sai_bi_ - c(l -bz+ Uz). 



Remembering that b± = ^\ — b\, we find ^[(AS)q — 

c{AB)i\ = —sa±b^^ < 0, and therefore that 
■^[{AB)q — c{AB)l] decreases with b^- Now, observ- 
ing that -^[{AB)q - c{AB)L]\b,=a,{i) = 0, we con- 
clude that for 5;, > 0^(1), ■^[{AB)q - c{AB)l] < 
0. Hence, {AB)q — c{AB)l decreases monotonically 
for bz G [02(1),!], and therefore takes values between 
[{AB)q - c{AB)L]\b^=i = (1 - c)az > -{I - c) and 
[{AB)q - c{AB)L]\b,=a,{i} = 1 - c. It follows that 
\{AB)q - c{AB)l \ < 1 - c, i.e. that (dJ) is satisfied. 

The case bz < az{—l) can be analyzed in a very sim- 
ilar way. We thus conclude that when Alice and Bob's 
settings are such that x = 0, but with b ^ Ha, the BGS 
model gives a valid EPR2 decomposition for all states. 



Appendix E: Derivation of Pa{t) 

We consider the case x = Oi assume for now that 
b lies in the Hardy sector Ha of a. The correlator reads 



i^B)l=^\ = J dtpaim - - 5,1). 



(El) 



To determine Pait) we will demand that condition (jlip 
is exactly saturated. 



Proposition 1. // Pa{t) 



■^'"(-^^ , then 

4c ^1-G(t)2 ' 



Proof. We will use the ansatz, for i G [— 1, 1], 

Pait) = [l + G{t)az]pit) (E2) 

where p{t) is a (normalized) even probability distribu- 
tion, and we further assume that 



±1 



dtp(t)Git) = l^ (4^(1; 



(E3) 



The ansatz (|E2[) ensures that pa (t) is a well-defined prob- 
ability distribution and that (A) ^ = a, as desired (which 
follow from the facts that 1 -I- G{t)az > and that 
p{t)G{t) is an odd function), while (|E3|) ensures that the 
proposition is true for a = b. At the end, we will have to 
check that condition (|E3[) holds and that p{t) is indeed a 
probability distribution. 



as 



We have 



bz = az{T) for some T £ [-1, 1]. 



i 

(^^)ir-f>=l-/ dtpa{t)\a,it)~bz\ 



(E4) 



(E5) 



=i+y" dtf{t)- j dtf{t) (E6) 

with f{t) = Pa{t)[az{t)-b,] 

p{t)[az — bz + {1 — azb2)G{t)]. By decomposing the 

two integrals, we find 



iAB)l7°^=l+ dtf{t) + 2 dtfit)- dtf{t) 
J -I Jo Jo 

=1- f dt[f{t)- f{^t)]+2 r dtf{t) 

Jo Jo 

=l-{l-azbz)^—^+2 [ dtfit), (E7) 
c Jo 

where we made use of the parity of p{t) and of assump- 
tion (lE3t. 



Demanding that {AB)^'" - c{AB)\-^ = 1 - c with 



{AB)^-° =azb,,+ sax_bi 



(E8) 



leads to 



2c 



: j dtfit) = sa^b^ - s(l - azb,). (E9) 


Expressing now bz as j^:^fja~^ S^^^ 



1 + G(T)a^ 



1 - a^b. 



1 -f- G(T)a, 
2 

1 + G(T)a. 



-G(r), 



(ElO) 



Using these expressions, (jE9[) leads to the integral equa- 
tion 

dip(t)[G(T) - Git)] = ^ (1 - Vl-G(T)2) , 

(Ell) 

which can be solved easily: e.g., after substituting 

t Git), (E12) 
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one gets an equation with a linear kernel, 
solution reads: 



s ln(7) 



1 



4c 



The unique 



(E13) 



As can be checked, p{t) is an even, normalized probability 
distribution on the interval [—1,1], and condition (|E3[) 
holds. 

□ 



The choice pa{t) 



sln(7) 1+G(t)a^ 



. , thus ensures that 

our maximal EPR2 decomposition, with pL = PLC = c, 
is valid for all measurement directions a and b such that 
X = 0, when b lies in the Hardy sector Ha of a. If x = 
but b ^ Ha, our local model gives the same correlation as 
the BGS model, which we proved in Appendix ID] to also 
give a valid EPR2 decomposition. Note that if x = tt, 
the situation is quite similar to x = 0, and one can prove 
again that our EPR2 decomposition is still valid. 

The question now is whether our maximal EPR2 de- 
composition remains valid for all other settings, when 
< Ixl < TT. The difficulty in checking that ([TT|) indeed 
holds is that one needs to calculate (^i?)L|a,b- given by 



{AB) 



L\a,b 



dtpa{t)EE''''M),Kx), (E14) 



with 



E 



2 / ax&2 - a^6_LCOsx 

—Uz arctan j — j 

TT V 6^sm|x| 



a±bz cosx 



2 , ^ f azb_i_ 

—Oz arctan 

TT \ a_LSm|x 



(E15) 



from the BGS model [28|. 

Unfortunately we were not able to calculate the in- 
tegral (jE14p explicitly. However, we carried intensive 
numerical checks to convince ourselves that ([TT]) always 
holds (up to machine precision of absolute order 10~^°), 
and therefore that our maximal EPR2 decomposition is 
valid for all possible measurement settings. 



Appendix F: Maximal violation of the CHSH 
inequality of the nonlocal part 



The correlator of the nonlocal part Pn l of our EPR2 
decomposition is given by: 



ENL{a,b) := (^B)^^, 



1 - c 



We denote with 5max the maximum violation of the 
CHSH inequality 



max ENL{ai,bi) + ENL{ai,b2) 

+ENL[a2M)- ENL{a2M)[ (Fl) 

/3max 




FIG. 3. Maximum violation /3max of the CHSH inequality by 
the nonlocal part Pnl of our EPR2 decomposition. 



In figure [3l the numerically determined values of /3max 
arc plotted for all states 6 G [0, 7r/2]. All states with 
6 < ^ violate Tsirelson's bound 2^2 [52|, approaching 
the algebraic maximum of 4 for — > 0; this proves that 
the corresponding nonlocal parts Pml arc nonquantum. 



